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Quantum theory predicts the existence of the 
Casimir force between macroscopic bodies, due 
to the zero-point energy of electromagnetic field 
modes around them. This quantum fluctuation- 
induced force has been experimentally observed 
for metallic and semiconducting bodies, although 
the measurements to date have been unable to 
clearly settle the question of the correct low- 
frequency form of the dielectric constant disper- 
sion (the Drude model or the plasma model) to 
be used for calculating the Casimir forces. At 
finite temperature a thermal Casimir force, due 
to thermal, rather than quantum, fluctuations of 
the electromagnetic field, has been theoretically 
predicted long ago. Here we report the experi- 
mental observation of the thermal Casimir force 
between two gold plates. We measured the at- 
tractive force between a flat and a spherical plate 
for separations between 0.7 /^m and 7 /im. An 
electrostatic force caused by potential patches on 
the plates' surfaces is included in the analysis. 
The experimental results are in excellent agree- 
ment (reduced % 2 of 1.04) with the Casimir force 
calculated using the Drude model, including the 
T = 300 K thermal force, which dominates over 
the quantum fluctuation-induced force at separa- 
tions greater than 3 /im. The plasma model result 
is excluded in the measured separation range. 

There are four known fundamental forces: electro- 
magnetism, gravity, weak, and strong interactions. The 
weak and strong interactions manifest themselves on 
length scales on the order of the size of a nucleus, at 
larger distances electromagnetism and gravity prevail. 
It may therefore come as a surprise that two macro- 
scopic non-magnetic bodies with no net electric charge 
(or charge moments) can experience an attractive force 
much stronger than gravity. This force was predicted 
by Hendrik Casimir in the late 1940s, and now bears his 
namei. The existence of this force is one of the few direct 
macroscopic manifestations of quantum mechanics, oth- 
ers are superfluidity, superconductivity, kaon oscillations, 
and the black body radiation spectrum. 



The first experimental test of Casimir's prediction 
came within a few years^, but the first "precision" mea- 
surements were made in the late 1990s^, leading to re- 
newed interest in this subject, both for experiment and 
theory 4 -^. Casimir's result, which he obtained by adding 
up the zero-point energies of the electromagnetic field 
modes of the cavity comprising two perfectly-conducting 
flat plates of area A separated by distance d, gives the 
force between the two plates at zero temperature: 
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Note the presence of the Planck's constant h, which in- 
dicates that this is a quantum effect, which vanishes 
in the classical limit. With some notable exceptions 6 , 
most Casimir force experiments, including ours, are per- 
formed with one flat plate and one spherical plate, to 
avoid the difficulty of aligning two flat plates to be par- 
allel to within fractions of a microradian. In the sphere- 
plane configuration, the zero-temperature Casimir force 
for perfect conductors can be deduced from equation (|TJ) 
using the proximity force approximation (valid when 
d < R)2A: 
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where R is the radius of curvature of the spherical plate, 
and d is defined as the distance between the flat plate 
and the closest point on the spherical plate. 

Real-world plates, however, are never perfectly con- 
ducting, but instead are characterized by a complex per- 
mittivity e(w) = e'{oS) + ie"(u>), a material-dependent 
function of frequency uj. This modifies the zero- 
temperature Casimir force, which now has to be cal- 
culated in the framework of the Lifshitz theory^, based 
on the computation of the electromagnetic field stress 
tensor, taking into account the correlated fluctuating 
charges and currents in the plates. The Lifshitz theory 
has been tested in experiments with < 200 A-thick liquid- 
helium films on cleaved crystal surfaces^, and measure- 
ments of the short-range forces between plates of a num- 
ber of metallic 11-13 and semiconducting 14 materials have 
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been found to be in reasonable agreement with the Lif- 
shitz theory at plate separations less than 1 ^m. Repul- 
sive forces have been observed between metallic plates 
immersed in fluids, also in agreement with the Lifshitz 
theoryi£i£. 

When the Casimir force between a pair of plates at 
finite temperature is considered, another source of the 
force arises - thermal fluctuations of the electromagnetic 
field. According to Bose-Einstein statistics, the popula- 
tion of an electromagnetic field mode of frequency u> at 
temperature T is 
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where the first term describes the zero-point fluctuations, 
and the second term is the thermal population for a 
gas of Bose particles (photons); k B is Boltzmann's con- 
stant. This second term gives rise to a finite-temperature 



term in the Casimir force F. 
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At small separations 



(rf < 1 /im at T = 300 K) the thermal force is much 
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smaller than the zero-point force: F<\ ; -C F^, , but 
for large separations (d > 3 /j,m at T = 300 K) the ther- 



mal force dominates: Fq ^> F^f ~"' . Although an anal- 
ogous thermal Casimir-Polder force has been measured 
between an atom and a surface^, the thermal Casimir 
force between two macroscopic objects has not been ex- 
perimentally observed until now. 

(T) 

The magnitude of the force F^, has been a subject of 
intense theoretical debate^—. The controversy concerns 
the form of the behaviour of the plates' complex permit- 
tivity e(ui) at low frequencies, and the related question of 
whether or not the zero-frequency transverse-electrical 
(TE ui = 0) mode should be included in the calcula- 
tion of the finite-temperature force. If the Drude model 
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frequency and 7 is the dissipation rate) is used to de- 
scribe the low-frequency permittivity, then the TE uj = 
mode does not contribute to the finite-temperature force, 
whose magnitude at large plate separations is calculated 
to be (for a sphere-plane geometry with d -C R) 
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where £ is the Riemann zeta function, £(3) ~ 1.202. The 
use of the Drude model, however, has been claimed to 
contradict the third law of thermodynamics (Nernst's 
heat theorem), as it appears to give rise to finite zero- 
temperature entropy of the electromagnetic field between 
the plates, if the plate material dissipation is allowed 
to vanish in the limit T — > 0, as is the case, for ex- 
ample, in an ideal crysta l 19 ' 20 (if the dissipation stays 
finite as T — > 0, there is no contradiction with the third 
law). In order to avoid this problem, a different form 
of the plates' low-frequency permittivity has been pro- 
posed: e p iasma(w) = 1 — Wp/w 2 , this is known as the 
plasma model. In this model, the TE uj = mode does 



contribute to the finite-temperature force, whose magni- 
tude at large plate separations is double the Drude result: 
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Some Casimir force measurements at small plate sep- 
arations (d < 750 nm) have been interpreted as be- 
ing in agreement with the plasma model 2 ^, although at 
such small separations the relative difference between the 
Drude and the plasma model predictions is very small. 
A calculation of the entropy of the electromagnetic field 
in the plasma model gives zero at absolute zero temper- 
ature, in agreement with the third law. Let us note, 
however, that the third law of thermodynamics does al- 
low for a non-vanishing entropy at zero temperature for 
systems with a degenerate ground stated, as is the case 
with glasses 2 ^. It is possible that the Drude model with 
vanishing dissipation in the limit T — > leads to a degen- 
erate ground state of the system due to persistent eddy 
currents (Foucault glass) 2 ^. Equations ((4]) and §5§ dis- 
play only the leading temperature-dependent term in the 
Casimir force; at 300 K and at the larger plate separa- 
tions considered in our experiment this is the dominant 
term. 

We report the measurement of the finite-temperature 
contribution to the Casimir force between gold plates. 
The total force between the plates measured in our ex- 
periment can be written as the sum of the Casimir force 
Fc(d) and the electrostatic force: 
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where eo is the permittivity of free space, and V is 
the computer-controlled bias voltage applied between the 
plates. The "minimizing potential" offset V m is due to 
the contact potential difference of approximately 20 mV 
between the two plates, caused by the several solder 
contacts around the electrical loop connecting the two 
plates. The second term in brackets is due to the re- 
gions (patches) of varying potential on the plate surfaces, 
caused, for example, by spatial changes in surface crys- 
talline structure, adsorbed impurities, or oxides, and in- 
variably present even on chemically inert metal surfaces 
prepared in an ultra-clean environmen t 25 ' 26 . The force 
caused by these patches has been experimentally ob- 
served in, for example, Refsi 14 ' 27 , and is characterized by 
the parameter V rms , related to the magnitude of the volt- 
age fluctuations across the plates. There are three length 
scales relevant to the form of the electrostatic patch force: 
the plate separation d, the "effective interaction length" 
r e ff = V Rd, and the typical patch size A. As shown 
in Refj2£, potential patches of size A <C d lead to an 
exponentially-suppressed electrostatic force between the 
plates (oc e~ d / x ), which we neglect. Potential patches of 
size A > r e f t give rise to an electrostatic force of the form 
ire R(V m (d) + Vi) 2 /d, where V m (d) describes the depen- 
dence of the minimizing potential on plate separation d, 
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and V\ is a constant^?.. We find that for the gold-coated 
plates used in our experiment V m is very nearly indepen- 
dent of d (variation is 0.2 mV between 0.7 /xm and 7 /im), 
therefore this additional force is small compared to the 
experimental error, and we do not include it in our data 
analysis. Finally, potential patches on the length scale 
d<A< r e f f give rise to an electrostatic force between 
the plates of the form V 2 ms /d, which is the last term in 
Eq. ©. 
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FIG. 1: The top-view schematic of the torsion-pendulum ex- 
perimental apparatus. 

A top-view schematic of our torsion pendulum appa- 
ratus is shown in Figure [1] and described in detail in the 
Methods section. The total force between the two plates 
is measured at 30 logarithmically-spaced plate separa- 
tions between 0.7 /im and 7 /^m in a series of 383 sweeps, 
adding up to a total of 8 days of data taking. In order to 
find the plate separation d, the force as a function of bias 
voltage V is recorded at a fixed d. A parabolic fit to the 
force- vs- voltage data is used to extract the separation d, 
and the offset V m - This procedure is repeated at sepa- 
rations of 0.7 /im and 7 /im, at intermediate separations 
the bias voltage is set to V m , eliminating the first term 
in brackets in equation (|6]), and d is determined from 
the change in the piezoelectric transducer strain gauge 
reading, pre-calibrated in a separate series of direct mea- 
surements. The closest approach of 0.7 /im was set due to 
feedback instability at smaller plate separations, caused 
by the large force gradient 3 -. 

A systematic correction has to be applied to the data 
to take into account fluctuations in plate separation d 3 °. 
The sources of these fluctuations are surface roughness of 
the plates, and pendulum fluctuations, caused, for exam- 
ple, by apparatus vibrations. Surface roughness measure- 
ments were performed with the Micromap TM-570 inter- 
fcromctric microscope at the Advanced Light Source Op- 
tical Metrology Laborator y 31 ' 32 , yielding an rms rough- 
ness of S q « 10 nm for the curved plate, and S q « 1 nm 
for the flat plate. Vibration-caused fluctuations in d were 
measured by connecting an inductor in parallel with the 
Casimir plates, and monitoring the resonance frequency 
of the resulting LC-circuit; rms fluctuations of < 40 nm 



were recorded. In addition, statistical error of ±10 nm in 
determination of d contributes in quadrature to the fluc- 
tuations mentioned above. We take the total rms plate 
separation fluctuation of S — (40±20) nm. From the Tay- 
lor expansion of the Casimir force about the mean plate 
separation, we deduce that a correction term F^S 2 /2 
has to be added to the theoretical force when compar- 
ing with experiment, the double prime denotes second- 
order derivative with respect to d. In addition, since the 
same correction exists for the electrostatic force, the plate 
separation d extracted from the electrostatic calibration 
was corrected by a factor 1 + (5/d) 2 , and the electrostatic 
patch force Vj? ma /d was corrected by the same factor. 
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FIG. 2: Experimental results for the total short-range force 
between gold plates. The data have been binned for clarity, 
the vertical error bars include contributions from the statis- 
tical scatter of the points as well as from uncertainties in 
the applied corrections, discussed in the text. Also shown 
are the four theoretical models for the Casimir force: (i) the 
Drude model including the T = 300 K finite-temperature 
force (red), (ii) the plasma model including the T = 300 K 
finite-temperature force (green), (iii) the Drude model with- 
out the finite-temperature force, ie with T — (blue), and (iv) 
the plasma model at T — (magenta). The data are plot- 
ted with F x d on the y-axis, so that the electrostatic force, 
proportional to 1/d, appears as a constant offset on the plot. 
Inset: experimental data with each of the theoretical Casimir 
force models subtracted, color as above. Electrostatic patch 
force neoRVrms/d corresponds to a constant offset (up to the 
small l + (S/d) 2 correction). The fit to the Drude model points 
is shown by the black line. 

The binned force data, corrected as described above, 
are shown in Figure [2j These data were taken with the 
bias potential set equal to the offset V m (determined at 
7 /im), therefore, as shown in equation ([5]), the recorded 
force is a sum of the Casimir force and the electrostatic 
patch-potential force, given by the second term in the 
brackets. Unfortunately an independent measurement 
of this electrostatic force with the required accuracy is 
currently not feasible. In order to perform such a mea- 
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surement, the local surface potential has to be measured 
with millivolt sensitivity and micron spatial resolution. 
Commercial Kelvin probes (used for measurements in 
Refi2£, for example) lack the necessary spatial resolu- 
tion, since, in order to get an appreciable sample-tip 
capacitance, a tip size of at least 1 mm is needed 33 , 
and state-of-the-art custom-built Kelvin probes have 
achieved micron-scale resolution, but have the electro- 
static potential sensitivity of only 30 mV— . Therefore 
the electrostatic patch-potential force has to be modeled 
and extracted from our data. This is done by fitting 
the experimental data with the expression of the form 
F(d) = F c (d) + ne RV 2 ms /d + a, where F c {d) is the 
theoretical prediction for the Casimir force with no ad- 
justable parameters (see below), and the constant force 
offset a is due to voltage offsets in measurement elec- 
tronics; for clarity we subtracted this offset from the dis- 
played data. The two fit parameters V rms and a are the 
only adjustable parameters used in our data analysis. 

We consider four theoretical possibilities for the 
Casimir force between the gold plates: (i) the Drude 
model including the T — 300 K finite-temperature force, 
(ii) the plasma model including the T = 300 K finite- 
temperature force, (iii) the Drude model without the 
finite-temperature force (ie with T = 0), and (iv) the 
plasma model at T = 0. The Casimir force for each of 
these models is calculated using the Lifshitz formalism 
and the gold optical permittivity data^, extrapolated to 
zero frequency using the corresponding (Drude/plasma) 
model with the parameters uj p = 7.54 eV, 7 = 0.051 eV 
(details in Methods section). The resulting theoretical 
force curves are shown in Figure^] Note that the data are 
plotted with F x d on the y-axis, so the patch-potential 
force contribution leads to a separation-independent off- 
set (up to the small 1 + (S/d) 2 correction) of the experi- 
mental data points compared to the theoretical models. 

The fit F(d) = F c {d) + ire RV 2 ms /d + a is performed 
independently for the four theoretical models. The best 
agreement is obtained for the Drude model at T — 300 K, 
with the reduced X 2 of 1.04. The rms patch potential 
fluctuation obtained from the fit is V rms = (5.4±0.1) mV, 
consistent with the magnitude of potential fluctuations 
expected across a gold surface due, for example, to work 
function variations^. The force offset obtained from the 
fit is a = (—3.0 ± 0.4) pN. Having extracted the elec- 
trostatic contribution to the force, and the offset, we can 
subtract them from the experimental data, the remainder 
is the Casimir force, plotted in Figure [3] together with the 
theoretical prediction of the Drude model at T = 300 K. 
We plot F x d 2 on the y-axis, so that the 1/d 2 finite- 
temperature force of equation ((4]) corresponds to an off- 
set that dominates the force at large plate separations 
(d ^ 3 /im). The experiment is in excellent agreement 
with the Drude model containing the thermal Casimir 
force. Note that the theoretical calculation of the ther- 
mal Casimir force, that dominates at large separations, 
is largely independent of the exact values of the Drude 
parameters. 
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FIG. 3: The short-range force data corrected for an elec- 
trostatic force with V rm s = 5.4 mV. The reduced x 2 °f 
1.04 demonstrates excellent agreement with the Drude model 
including the thermal Casimir force at T — 300 K (red 
lines). The grey band represents theoretical uncertainty in the 
Casimir force calculation from the ellipsometry data; the force 
curve with Drude parameters u p = 7.54 eV, 7 = 0.051 eV was 
chosen for comparison with experiment (see Methods section). 
The data are plotted with F x d 2 on the y-axis, so that the 
thermal Casimir force, given by equation @, corresponds to 
an offset of 97 pNx/im 2 , which dominates the force at large 
plate separations. In this region the Casimir force is largely 
independent of the material properties of the plates. 

The experimental data rules out the other three the- 
oretical Casimir force models. The reduced % 2 obtained 
from the fit is: (ii) \ 2 = 32 for the plasma model at 
T = 300 K (best fit V rms = 3.0 mV), (iii) X 2 = 23 for 
the Drude model without the T = 300 K temperature 
correction (best fit V rms — 4.0 mV), and (iv) x 2 = 43 
for the plasma model without the T = 300 K tempera- 
ture correction (best fit V rms = 3.6 mV). Therefore our 
experiment rules out the plasma model for the Casimir 
force between gold plates in the separation range 0.7 /im 
to 7 /im, confirms the Drude model, and demonstrates 
the existence of the T = 300 K thermal Casimir force. 
The thermal Casimir force drops off as T/d 2 , and there- 
fore dominates over the purely quantum T = Casimir 
force (which behaves roughly as 1/d 3 ) for plate separa- 
tions greater than 3 /im. 



I. METHODS 

A torsion pendulum is suspended inside a vacuum 
chamber (pressure 5 x 10 -7 torr) by a tungsten wire of 
25 /im diameter and 2.5 cm length. The force to be mea- 
sured is between the two "Casimir plates", each coated 
with a 700 A (optically thick) layer of gold evaporated 
on top of a 100 A-thick layer of titanium. One is a flat 
plate mounted on one side of the pendulum, as shown in 
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Figure [TJ the other is a spherical lens (radius of curvature 
15.6 cm, measured with the Micromap TM-570 interfero- 
metric microscope at the Advanced Light Source Optical 
Metrology Laboratory together with surface roughness, 
and found to vary by less than 2% over the surface of 
the lens), mounted on a Thorlabs T25 XYZ positioning 
stage, which, together with a piezoelectric transducer, is 
used to vary the plate separation d. The attractive force 
between the plates creates a torque on the pendulum 
body, which is counteracted by a pair of "compensator" 
electrodes on the opposite end of the pendulum. The 
voltage that has to be applied to the compensator elec- 
trodes to keep the pendulum stationary is proportional 
to the force between the Casimir plates, with the calibra- 
tion coefficient extracted from the measurements of the 
electrostatic force between the plates. Pendulum rota- 
tion is detected by a capacitance bridge connected to the 
compensator electrodes, and a servo loop is used to ap- 
ply the compensator electrode voltage necessary to hold 
the pendulum in equilibrium (at a fixed angle). A Nd- 
FeB magnet is placed under the pendulum body to damp 
the swinging modes of the pendulum. The experiment is 
placed on a vibration-isolation slab, extending down to 
the bedrock below the building foundation. Subsequent 
measurements indicated that a smooth long-term drift, 
which is subtracted from the experimental data, is due to 
tilting of this slab (measured by an Applied Geomechan- 
ics model 701-2 A tilt meter), correlated with ambient hu- 
midity. Both humidity and temperature were monitored 
using an Onset HOBO U12 data logger. Temperature 



variations were less than 1 C. 

Ellipsometric measurements in the wavelength range 
191 nm to 1700 nm were performed on gold-coated plates 
prepared identically to the ones used in the force mea- 
surement apparatus. The resulting absorption and reflec- 
tivity data were within a few percent of the data given 
in Ref. 35 . In order to calculate the theoretical Casimir 
force with the Lifshitz formalism, the data were extrapo- 
lated to higher frequencies using the data in ReL^ 5 ., and 
to zero frequency using either the Drude or the plasma 
model. The following ranges for the plasma frequency 
and dissipation were explored: lj p = 6.85 eV to 9.00 eV, 
and 7 = 0.02 eV to 0.061 eV, resulting in theoretical 
force variation of up to 3%, shown by the grey band in 
Figure [3l The parameters uj p — 7.54 eV, 7 = 0.051 eV 
were chosen for the theoretical Casimir force that was 
compared with experimental data, since they best fit the 
optical data given in Ref.—. 
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